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ABSTRACT 

We present observations of Neptune's 1- and 3- mm spectrum from the Combined Array for Research 
in Millimeter-wave Astronomy (CARMA). Radiative transfer analysis of the CO (2-1) and (1-0) 
rotation lines was performed to constrain the CO vertical abundance profile. We find that the data are 
well matched by a CO mole fraction of O.I^q' 1 parts per million (ppm) in the troposphere, and I.I^q' 2 . 
ppm in the stratos phere. A flux of 0.5- 20 x 10 s CO molecules cm -2 s _1 to the upper stratosphere is 
implied. Using the Za hnle et al.l (|2003l ) estimate for cometary impact rates at Neptune, we calculate 
the CO flux that could be formed from (sub)kilometer-sized comets; we find that if the diffusion 
rate near the tropopause is small (200 cm 2 s -1 ), these impacts could produce a flux as high as 
0.51q'4 x 10 s CO molecules cm -2 s -1 . We also revisit the calculation of Neptune's internal CO 
contribution using revised calculations for the CO— >CH4 conversion timescale in the deep atmosphere 
(|Visscher and M oses 20Q]). We find that an upwelled CO mole fraction of 0.1 ppm implies a global 
O/H enrichment of at least 400, and likely more than 650, times the protosolar value. 

NOTICE: this is the authors' version of a work that was accepted for publication in Icarus. Changes resulting 
from the publishing process may not be reflected in this document. A definitive version was subsequently 
published in Icarus 222 (2013) 379-400. 



1. INTRODUCTION 

In equilibrium, carbon monoxide (CO) should be con- 
fined to the warm interiors of the Solar System giant 
planets. Therefore, its detection in the upper atmo- 
spheres of all four of these planets (lBeerlll975t Noll et all 



ll986HEncrenaz et al.ll2004HMarten et a 1.11 199 11) indicates 
disequilibrium processes at work. Two pathways exist for 
enriching the atmosphere in CO: vertical mixing from the 
deep atmosphere and external supply from the environ- 
ment. 

CO production occurs via the net thermochemical re- 
action 



CH 4 + H 2 = CO + 3H 2 . 



(1) 



At the temperatures and pressures of Neptune's atmo- 
sphere, the left-hand side of this equation dominates, 
with nearly all of the carbon present in the form of 
CH4. Indeed, CH4 has been observ ed in Neptune's tro - 
posphere at an abundance of 2.2% (|Baines et all . 1995). 
CO is more stable and therefore more abundant at the 
warmer temperatures of the deep atmosphere. Conse- 
quently, CO abundances in the upper atmosphere can ex- 
ceed their equilibrium value if convective transport from 
the warm i nterior is more rap i d tha n the CO destruc- 
tion rate ()Prinn and Barshavl Il977t iFeglev and Prinnl 
fl986h . The CO abundance originating from the in- 
terior therefore depends on the internal water abun- 
dance (Eq. (1)) as well as the speed of vertical mixing, 
and there by acts as a chemical probe of the deep at- 
mosphere (IFeglev and Loddersll9 94; Lod ders and Feglevl 
[1991 [2002TlVisscher et al.ll2010l) . 

An alternative source of CO in giant planet atmo- 

1 Astronomy Department, University of California, Berkeley, 
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spheres is the planetary environm ent. Observations 
of st ratospheric water and CO9 (IFeuchtgruber et al.l 
Il997t Ide Graauw et al.lll997t ILellouch et all 119971 ) indL 
cate an external supply of oxygen to the giant plan- 
ets, which could originate from rin gs and icy satellites 
or interplanetary dust (me teoroids) (|Feuchtgruber et all 
Il997t lEncrenaz et al.lU999h . Oxygen then forms CO by 
combining in the st ratosphere with the byproducts of 
meth ane photolysis (jRosenqvist et al.11199 2; Mose s et al.l 
2000). However, the CO/H 2 ratios observed for Jupiter 
and Neptune seem to be inconsistent with oxygen supply 
by these mechanisms (ILellouch et al.l 2002; Bczar d~et al.l 
I2002t ILellouch et all 120051). Shock chemistry due to in- 
fall of comets also can produce significant amounts of 
CO in the upper atmosphere, as demonstrated by the 
1994 impact of come t Shoemaker-Levy 9 with Jupiter 
(jLellouch et al.l I1995D . Infall of (sub)kilometer-sized 
comets has been shown to be sufficient for supplying 
Jupiter's observed stratospheric CO abundance if the 
eddy mixing coeffi cient near the trop opause is smaller 
than 300 cm 2 s" 1 (jBezard et all 120021) . Once it is pro- 
duced, CO is stable and is removed from the stratosphere 
by downward transport. 

Atmospheric CO enrichment by these two very differ- 
ent pathways can be distinguished by measuring the ver- 
tical CO profile: a uniform distribution of CO through- 
out the upper atmosphere indicates that CO is be- 
ing mixed up to observable levels from the deep at- 
mosphere. If instead CO is being produced in the 
stratosphere of the planet, then downward transport will 
act as a sink, and there will be a higher CO abun- 
dance in the stratosphere than the tropo sphere. Early 
meas u rements of Neptune's stratospheric (IMarten et al.l 
119931: IRosenqvist et all Il992t IMarten et al.l 120051) "and 
tropospheric (|Guilloteau et al.lTl993l: INavlor et all 119941: 
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iCourtin et al]|1996t lEncrenaz et~aIlll996D CO mole frac- 
tions were all roughly consistent with a CO abundance 
of 1 part per million (ppm), which led several authors 
to tentatively conclude that Neptune's CO was interna l 
in origin (e.g. ICourtin et aDll996t [Marten et al.1 120051 ). 
However, the estimates of Neptune's CO abundance from 
these observa tions were highly divergent. More recent 
observations (ILellouch et all 120051: IHesman et~aT1 12007ft 
suggest that both internal and external sources play a 
role in supplying CO to Neptune: these experiments 
simultaneously determine the stratospheric and tropo- 
spheric abundances from the shapes of CO rotational 
lines. Figure 1 illustrates the contributions of the var- 
ious atmospheric levels to the line intensity at a range 
of frequency offsets from line center, for the three low- 
est CO rotational transitions. Peaks in the contribution 
functions range from above 0.1 mbar at line center, down 
to several bars in the far wings indicating that charac- 
terization of the full CO vertical profile requires both 
high frequency resolution to measure the shape of the 
narrow (~10 MHz wide) central line peak, and broad 
frequency coverage (of order 6-10 GHz) to probe the 
CO mole fraction below the tropopause. From their 
analys es of the CO ( 2 -1) a nd C O (3-2) lines, r espec - 
tively, iLellouch et all (|2005[ ) and IHesman et all (|2007f ) 
conclude that Neptune's CO mixing ratio is measur- 
ably higher in the stratosphere, indicating that CO has 
both an internal and external origin. Such a dual origin 
is also indicated for Jup iter and Saturn ([Bezard et al.1 



2002t iCavalie et aTll2009fl . Quantitatively , however, the 
Lellouch et al.l (|2005[) and IHesman et all (|2007f ) Neptu- 



nian CO profiles are inconsistent to within their quoted 
uncertainties, particularly in the upper stratosphere. Un- 
certainties in calibration and in Neptune's thermal profile 
may be responsible for discrepancies between these pub- 
lished CO values. Recent observations at infrared wave- 
length s have not resolved the discrepancy: Fletcher et al. 
(2010) found a CO mixing ratio of 2.5 x 10~ 6 (mole frac- 
tion of 2.1 ppm) a t altitudes above 10 mbar, which is 
consistent with the Hesma n et al.1 (|2007[ ) result, whereas 
ILellouch et al.1 (|201Q[ ) favo r a 1 ppm stratosphe ric mole 
fraction in agreement with ILellouch et all (|2005h . 

In this paper we present new observations and model- 
ing of the CO (2-1) and (1-0) rotational lines in Nep- 
tune's millimeter spectrum in order to better constrain 
the vertical CO profile. We investigate how the un- 
certainty in Neptune's atmospheric thermal profile af- 
fects these constraints. This is followed by a discus- 
sion of the implications of our derived CO profile on 
Neptune's global oxygen abundance: updated labora- 
tory measurements of reaction rate s have shown that 
the chemical scheme proposed bv iPrinn and Barshavl 
(|1977l ) and adopted bv lLodders and Feglevl (|1994T > is ac- 
tua lly too slow to be relev a nt for CO quenching kinet- 
ics (|Griffith and Yelldfl999h . IVisscher and MosesJ (|20H 
have further revised the kinetic scheme for CO— >CH4 
conversion using new values for the reaction rate coef- 
ficients. Using this new rate-limiting step, we calculate 
the CO mole fraction that is transported upwards from 
Neptune's deep atmosphere, ass uming an effe ctive mix- 
ing leng th scale as de t ermin ed bv lSmithl ((19981 ) . We also 
use the iZahnle et al.l (|2003l ) impact rates to determine 
the effectiveness of cometary impacts in supplying oxy- 
gen to Neptune's atmosphere for stratospheric CO pro- 
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Fig. 1— Contribution functions for the CO (1-0), (2-1), 
and (3-2) lines, illustrating the altitudes contributing to the 
line at offsets from to 4.5 GHz from line center. The contri- 
bution functions include H2-H2, HVHe and H2-CH4 collision- 
induced absorption, and absorption due to a constant 1 ppm 
abundance of CO. Close to line center, the emission originates 
in the stratosphere; at larger offsets from line center, deeper 
levels contribute to the line. The CO (1-0) and (2-1) lines 
are the subject o f this work; th e CO ( 3-2) line was observed 
most recently by Hesman " et al.1 l|2007f l. 
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duction. 

2. DATA 
2.1. Observations 

Disk-integrated observations of Neptune in the J = 1 
and J = 2-1 transitions of CO were performed with 
the Combined Array for Research in Millimeter- wave As- 
tronomy (CARMA) in March and April 2009. CARMA 
is a 23 element interferometer that combines six 10-meter 
antennas, nine 6-meter antennas, and eight 3.5-meter an- 
tennas. Our observations were performed with the 6- and 
10-meter antennas only, for a total of 105 baselines and 
2900 m 2 of total collecting area. CARMA rotates be- 
tween 5 different standard antenna configurations; our 
CO (2-1) data were taken in CARMA's D array which 
consists of baselines of 11-148 meters for a synthesized 
beam of 3.0x1.9 arcseconds at 230 GHz. Our CO (1-0) 
data were observed in CARMA's C array, with baselines 
of 26-370 meters for a synthesized beam of 2.1x1.8 arc- 
seconds at 115 GHz. In both sets of observations, the 
beam size is comparable to the size of the planet; there- 
fore we do not obtain information about spatial varia- 
tions in the CO distribution. 

At the time of our observations the CARMA correlator 
had three dual bands (or "windows") with configurable 
bandwidth of either 500, 62, 31, 8 or 2 MHz. Each band 
could be placed independently anywhere within the 4 
GHz IF bandwidth, and appears symmetrically in the 
upper and lower sidebands of the first local oscillator. 
The sideband, total bandwidth, and central frequency of 
the 6 windows (3 in each sideband) for each of our cor- 
relator setups are given in Table [l] In each setup, we 
tuned the receivers to the CO (1-0) or (2-1) rest fre- 
quency in either the lower (setup a,c) or upper (setup 
b,d,e) sideband. Band 2 was always configured with a 
bandwidth of 62 MHz across 63 channels, to observe the 
CO line center at a resolution of 0.98 MHz per channel. 
Bands 1 and 3 were each configured for maximum band- 
width (500 MHz across 15 channels) and positioned at 
a frequency offset from the CO line center. Bands 4-6 
correspond to Bands 1-3 in the opposite sideband. For 
the CO (1-0) line observations, we configured Band 1 at 
an offset of —0.75 GHz and Band 3 at an offset of +0.30 
GHz (correlator setup 'e' in Table [J). For the CO (2-1) 
line, we alternated between 4 different correlator setups 
of the 500 MHz bands (Table [U setups 'a-d') to extend 
our frequency coverage of the broad CO line wings. Fre- 
quency offsets of the 500 MHz bands from line center 
vary from 0.15 — 9.3 GHz. 

Each individual observation consisted of a 3.6 — 6.8 
hour track, with a 15 minute observation of a passband 
calibrator, followed by a series of observing cycles of 15 
minutes on source and 3 minutes on a phase calibra- 
tor. Optical pointing was performed every hour. No 
planets were available for primary flux calibration, so 
we observed MWC349 when it was available. The qual- 
ity of the data from tracks shorter than 3 hours and 
tracks with very poor weather were too low, and hence 
the data were rejected (2 of 8 tracks at both 1 mm and 3 
mm) . Weather conditions during the remaining observa- 
tions were generally fair to good, with typical root mean 
square (rms) path errors of 150-400 /mi on a 100-m base- 
line, and zenith optical depths of 0.1 0.2. The total 



time on source was about 14.9 hours at 1 mm and 18.6 
hours at 3 mm. These observations are summarized in 
Table H 

2.2. Calibration 

The data are reduced and calibrated using the 
MIRIAD software package (|Sault et all 1201 It ). Prior to 
any calibration, flagging is performed: we flag 10 edge 
channels for the narrow bands and 3 edge channels for 
the wide bands; the narrow window that is not centered 
on the emission core of the CO line (Band 5: see Ta- 
ble [I]); and any bad data. After performing passband 
calibration using our bright quasar observations, time- 
dependent gain solutions are derived using the wide-band 
data, and then applied to the full data set. We do an 
initial self calibration using the phase calibrator. This 
first calibration is performed in two steps: a record-by- 
record phase solution is found for the phase calibrator, 
to remove short-term variations. Then, a phase and am- 
plitude self calibration is performed using an 18 minute 
interval and applied to the Neptune data. Finally, a 
record-by-record phase-only self calibration is performed 
on Neptune itself to remove short-term phase variations 
in the Neptune data. 

2.2.1. CO (3=2-1) 

For this set of observations, only one track (28 March 
2009) contains observations of the primary flux calibra- 
tor, MWC349. We flux calibrate this track using an as- 
sumed value of 1.86 Jy for MWC349 at 230 GHz, based 
on flux density measurements obtained with CARMA be- 
tween 2007 and 20110. We estimate the flux calibration 
to be accurate to better than 20%. For a first pass at the 
flux calibration of the other tracks, we use our phase cal- 
ibrator as a secondary flux calibrator. We then bin the 
Neptune visibility data into 8 (u, v) bins between and 
80 kX, and adjust the flux of the phase calibrator to align 
the binned (u, v) data as much as possible: these adjust- 
ments are about 2-5%. From the deviation of these fits, 
we estimate the remaining error in day-to-day gains to 
be less than 3%. 

2.2.2. CO (3=1-0) 

All but one of our 3-mm tracks contains primary flux 
calibrator data (Table [2]). However, we find that in gen- 
eral, the flux determination from the primary flux cali- 
brator is very sensitive to bad/noisy data. We are able to 
get more stable fluxes by using our best primary flux cal- 
ibrator data from 28 April 2009, and applying the same 
procedure as for the 1-mm data. We use a flux density 
of 1.29 Jy for MWC349 at 3 mm 1 . 

2.3. Imaging 

Our 3-mm data sets are all observed with the same cor- 
relator setup; therefore, we combine all the 3-mm data 
prior to imaging. In contrast, the 1-mm tracks are ob- 
served at different frequencies and we image each track 
separately. After imaging each channel, we deconvolve 
the dirty maps using the CLEAN algorithm. Figure [D 
shows several of our final maps. Each map represents 

2 http: // cedarf lat .mmarray . org/f luxcal/pr imary_sp_index.htm 
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TABLE 1 

Correlator setups used in the observations. At 1 mm, each track uses one of setups 'a'-'d', while all 3-mm tracks use correlator 

setup 'e'. Bands ("windows") 1-3 are located in the same sideband as the CO line center frequency: Band 2 is always a 62 
MHz band centered at the CO line rest frequency, and Bands 1-3 are wide (500 MHz) bands offset from line center. Bands 4-6 
correspond to Bands 1-3 in the opposite sideband; these data are also recorded, and the wide-band data from both sidebands 

are used in our analysis. 



Configuration / setup 




a 


b 


c 


d 


e 


Band 1 


... . = 

sideband 
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— 1 run 

USB 
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Loo 


— I TO l S 

USB 


i — rrcT3 

USB 




oanawiatn ^IVlriz ) 


ouu 


ouu 


OUU 


OUU 


OUU 




center frequency (GHz) 


231.288 


230.318 


230.758 


229.538 


114.521 


Band 2^ 


sideband 


LSB 


USB 


LSB 


USB 


USB 




bandwidth (MHz) 


62 


62 


62 


62 


62 




center frequency (GHz) 


230.538 


230.538 


230.538 


230.538 


115.271 


Band 3 


sideband 


LSB 


USB 


LSB 


USB 


USB 




bandwidth (MHz) 


500 


500 


500 


500 


500 




center frequency (GHz) 


230.288 


229.218 


232.538 


230.688 


115.521 


Band 4 


sideband 


USB 


LSB 


USB 


LSB 


LSB 




bandwidth (MHz) 


500 


500 


500 


500 


500 




center frequency (GHz) 


235.788 


221.818 


239.258 


225.538 


110.021 


Band 5^ 


sideband 


USB 


LSB 


USB 


LSB 


LSB 




bandwidth (MHz) 


62 


62 


62 


62 


62 




center frequency (GHz) 


236.538 


221.598 


239.478 


224.538 


109.271 


Band 6 


sideband 


USB 


LSB 


USB 


LSB 


LSB 




bandwidth (MHz) 


500 


500 


500 


500 


500 




center frequency (GHz) 


236.788 


222.918 


237.478 


224.388 


109.021 



a ccntcred on the CO line 
b not used 



TABLE 2 

Summary of observations (all in 2009) 



Line 


Frequency (GHz) 


Date 


Correlator setupj^ 


T int (hours^ 


Calibrators 


rms path (fJ-va) 


Tzen 


CO (2-1) 


230.538 


07 Mar 


a 


2.33 


30454.3^,2229-085^ 


208 


0.22 






08 Mar 


b 


2.08 


3C454.3~, 2229-085^ 


165 


0.24 






10 Mar 


c 


2.83 


3C454.3 C , 2229-085 d 


157 


0.19 






11 Mar 


d 


2.33 


3C454.3 C , 2229-085 d 


133 


0.23 






27 Mar 


a 


2.57 


3C454.3 C , 2229-085 d 


162 


0.17 






28 Mar 


a 


2.71 


3C454.3 C , 2229-085 d , MWC349°_ 


147 


0.23 


CO (1-0) 


115.271 


16 Apr 


e 


4.25 


1751+096 c , 3C446 d , MWC349~ 


413 


0.12 






17 Apr 


e 


3.40 


3C446 c < d , MWC349 e 


306 


0.12 






23 Apr 


e 


3.03 


3C454.3 C , 3C446 d , MWC349 e 


335 


0.16 






24 Apr 


e 


2.25 


3C454.3 C , 3C446 d , MWC349 e 


424 


0.18 






28 Apr 


e 


3.15 


3C454.3 C , 3C446 d , MWC349 e 


200 


0.11 






29 Apr 


e 


2.50 


3C454.3 C , 3C446 d 


180 


0.10 



a as specified in Table 1 
b time on source 
c Passband calibrator 
d Phase calibrator 
c Flux calibrator 



a single channel of data: the top row are 1-mm chan- 
nel maps from a day with typical weather conditions (11 
March), and the bottom row are examples of 3-mm chan- 
nel maps. Below each map is a scan in Right Ascension 
through the center of Neptune. 

2.4. Flux determination and error estimate 

The flux density in each channel is determined as the 
integrated intensity over a hand-selected region, chosen 
to include all of the signal from Neptune, as indicated 
by the solid line in each map in Fig. [5] To quantify 
the errors in the determined flux density, we then shift 
this region to 8 different positions outside of the source 
(indicated by dashed lines in Fig. [5]). The error in each 
channel is estimated as the standard deviation of the inte- 
grated fluxes of these regions. We also calculate the flux 



errors using the rms of the pixel values in the residual 
maps. We find that the two methods of error estimation 
are in good agreement. These error estimates do not in- 
clude systematic effects, such as errors in the bandpass 
calibration or in the assumed flux of the calibrator; how- 
ever, they are useful weights for fitting the data to models 
(see Sections S] and O- Raw spectra, with uncertainties, 
are shown in Fig. [J2 for reference, the individual 1-mm 
tracks are differentiated by color. A frequency-binned 
version of the data is plotted in black, using a bin inter- 
val of three times the channel width. 

3. MODEL 

To model Neptune's millimeter spectrum, we devel- 
oped a line-by-line radiative transfer code that integrates 
the equation of radiative transfer, assuming local thermo- 
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Fig. 2. — Representative channel maps; color scale indicates intensity in Jy/beam. The top row are maps from a typical 1-mm 
day (11 March); the bottom row are the channel maps created from all days of 3- mm data. The left-most figures are wideband 
channels, roughly 5 GHz from line center. In the middle column are maps from wide channels in the absorption, about 500 
MHz from line center. Maps on the right are from narrow channels at line center. Below each map is a scan through the map 
Right Ascension (a value of on the y- axis). Neptune is unresolved in the maps; the total flux in each channel is determined 
by summing over the region indicated by the solid white line. The error is estimated by the regions indicated by dashed lines- 
see Section \2. 31 for discussion. 



dynamic equilibrium (LTE) 

B v {T D ,fj,)= / B„(T)e- T /»dT/fx (2) 
Jo 

over a model atmosphere consisting of 2000 plane-parallel 
layers extending from 200 bar to 5 /ibar. In this equa- 
tion, Tjj is the disk brightness temperature, /i = cos 9 
with 9 defined as the angle between the line of sight 
and local vertical, B V (T) is the Planck function for 
temperature T, and r is the optical depth. The code 
is optimized for least-squares fitting of the CO alti- 
tude profile (Section^]), and has been thoroughly tested 
a gainst the microwave r adiative transfer code described 
in lde Pater et al.1 (Il991af). after the latter was updated as 
described bv lde Pater etaTl (|2005t ). The Ide Pater et al.l 
(|1991al) code has been further revised to include the 



Orton et a l. (2007b) H2 absorption coefficients (Section 
3.21) . and absorption d ue to CO using subroutines from 
do P ater et all (|1991bD that were updated with the pa- 
rameters described in Section [3~2l We find that the two 
radiative transfer codes produce consistent model spectra 
for the same input parameters, except at the CO line cen- 
ter, because the new code includes the effect of Doppler 
broadening on the sh ape of the emi s sion pe ak (Section 
13. 4[) and the revised Ide Pater et al.l (|1991aD code does 
not. The details of the new model are described below. 

3.1. Composition 

N eptune's upper atmosp here is dominated by H2 and 
He. iConrath et al.l (| 19911) estimated the relative abun- 
dances of these two species using constraints on the at- 
mospheric mean molecular weight from Voyager infrared 
and radio occultation measurements; they found a best- 
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Fig. 3. — Raw CO (1-0) and (2-1) spectra. Each grey point corresponds to the flux density calculated from a single 3-mm 
channel map; colors are used to delineate between individual tracks at 1 mm. Error bars are calculated using the method 
described in Section 12. 31 Black points are binned to 3 times the width of the individual channels. 



fit helium mole fra ction of 0.19±0.032 . Th e detection of 
HCN on Neptune ()Marten et al.lll99lD led IMarten et al.l 
(119931) to suggest a scenario (originally proposed by 
iRomani et al lll989l ) in which Neptune's nitrogen is pre- 
dominantly present as N 2 , rather than NH 3 . Such a high 
N2 abundance is greater than expe cted from thermo- 
chemical equilibrium arguments (e.g. iFegley et al.lll99l"l 
for Uranus), but could be the source of at omic nitrogen 
for th e production of stratospheric HCN (|Marten et al.l 
119931) . and w ould help explain th e observed atmospheric 
NH 3 deficit (IRomani et al.l 119891 d iscuss ed below). Re- 
analy sis of the lConrath et ah HIT991 data (|Conrath et al.l 
1993) showed that a mole fraction of 0.003 of N2 and 0.15 
of He is consistent with the Voyager measurements. Us- 
ing spectra from the Infrared Space Observatory Long- 



Wavelength Spectrometer (ISO-LWS), iBurgdorf et al.1 
(pOOl found a He/(H 2 +He) mass ratio of 26.4+^%, and 
an N2 mixing ratio of less than 0.7%. They determined 
that the corr elation of their value of the He mole frac- 
tion with the iConrath et al.l (|1993l ) results gives an N 2 
mole fraction of 0.3 ± 0.2%. Accordingly, for our atmo- 
spheric models we maintain a He/H 2 ratio of 0.15/0.847 
and an N 2 /H 2 ratio of 0.003/0.847 by number through- 
out our model atmosphere; t hese numbers are co nsis- 
tent with the solution s of both IConrath et all ()1993l ) and 
IBurgdorf etaLl (|2003l) within the 1 a uncertainties. Since 
the presence of N 2 in Neptune's atmosphere is not cer- 
tain, we also tested models with the same He/H 2 ratio 
and no N 2 ; the difference was negligible. 
Previous observations have shown that CH4 is su- 
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persaturated in Neptune's stratosphere ()Orton et al.l 
Il987t IQrton et al.l [19901 lYelle et alJll99l . We adopt 
the stratospheric C H4 profile recently derived by 
iFletcher etaD (|2010T > from AKARI infrared data, which 
has a mole fraction of 9 x 10 -3 at 50 mbar. Below the 
methane condensatio n level, we adopt a m ole fraction of 
0.022 as measured bv lBaines et al.1 (|1995f l. which implies 
an en richment factor of ~ 50 over the protosolar C/H 
ratio (jAsplund et al.ll2009h . 

While the composition of Neptune's deep troposphere 
has yet to be uniquely determined, data from centimeter 
wavelengths suggest that NH3 is depleted in Neptune 's 
atmosphere (jRomani et al.lll989l : Ide Pater et al.l Il991af ) . 
Good fits to the cm wavelength range are obtained by 
enhancing H2S 30-50 times over the sol ar S/ H value 
(|de Pater et al.lll991aHDeboer and Steffeslll9961) and us- 
ing a s olar abundance of nitrogen in NH3 . Hoffman et al. 
(|2001f ) suggest that PH 3 may be an im portant source 
of mic rowave absorption as well; however, iMoreno et al.1 
(|2009f ) determined an upper limit of phosphorous in Nep- 
tune's upper atmosphere of 0.1 times the solar abun- 
dance, from observations of the PH3 (1-0) transition. 
Due to the uncertainty in our absolute flux calibration 
(as much as 20%, see Section [2]), we do not attempt to 
constrain the abundances of species other than CO from 
our data. However, we do investigate the potential effect 
of these absorbers on the 1- and 3-mm spectra: u sing 
the code described by Ide Pater et al.l (|1991aL 120051) , we 
model the deep atmosphere with a 10-50 times solar en- 
richment of H2O and H2S, a solar abundance of NH3, 
and 0.1 times solar PH3. Trace species are removed from 
the model atmosphere at higher altitudes by condensa- 
tion, when the partial pressures of the trace gases exceed 
the saturation vapor pressure (see Ide Pater et al.l (|2005h 
for a detailed description) . Figure 2] shows the effects 
on the microwave spectrum of absorption due to each of 
these gases individually, as well as that from all gases 
combined. 

3.2. Opacity 

To determine the optical depth t v (z) we consider the 
following sources of opacity: 

3.2.1. Collision-induced H2 absorption 

For Neptune, the dominant millimeter-wave opacity 
source is collision-induced absorption by H2 with H2 , He, 
and CH4. We use the absorpti on coefficients calcul ated 
from revised ab-initio models of lOrton et al.l (|2007b| ). as- 
suming an equilibrium distribution of hydrogen. These 
authors incorporate a correction to the Bq r ysow model s 
(|Borvsow et al.l[l985l [l988t IBorvsowlll99l Il992l U99l . 
and show that the new coefficients are an improvement 
at low temperature s. In practice, we fi nd that spectra 
produced using the IQrton et al.l ()2007bl ) coefficients are 
1.6-2 K higher than those made using the lBorvsow et all 
(1985) models. 

3.2.2. CO 

The absorption due to the CO (1-0) and (2-1) rotation 
lines is calculated assuming a Voigt line shape profile. 
The H2 broa dened line half- width is determined from a fit 
to the data of lMengel et~aTI (j2000D and is -2.8 MHz/Torr 
at 300 K for the CO (1-0) and (2-1) lines. The remain- 
ing line parameters are taken from the HITRAN 2008 



database (iRothman et al .1120091) . A Van Vleck-Weisskopf 
line sha pe was tested using the updated Ide Pater et al.l 
(|1991al) code and found to be a nearly identical match to 
the Voigt line shape in the atmospheric region probed. 

3.2.3. Additional microwave absorbers 

We model several other potential sources of radio- 
wavel ength opacity (H2O, NH 3, H 2 S, PH 3 ) using the up- 
dated Ide Pater et al.1 ( 1991aD code and the abundances 
from Section l3~Tl We find that H 2 and NH 3 do not af- 
fect the spectrum at 1-3 mm wavelengths (Fig. SJ. The 
far wings of the PH3 (1-0) line at 266.9 GHz could in- 
fluence the high frequency wing of the CO (2-1) line; 
however, for a 0.1 times solar abundance of PH 3 the line 
is too weak to significantly affect the CO profile. We 
note that a 0.1 times solar ab undance of PH3 is an up- 
per limit ([Moreno et al.ll2009h : larger values would also 
be inconsistent with the available data according to our 
own models, as shown. We find that absorption due to 
H2S is important at wavelengths longer than 1 mm; at 3 
mm the addition of > 10 times solar abundance of H2S 
decreases the continuum level by ~20 K. 

We therefore incorporate H2S absorption into our new 
ra diative transfer code follow ing the formalism described 
by iDeboer and Steffesl (119941); line paramete rs are taken 
from the JPL catalog ()Pickett et al.l 119981 ). while pa- 
rameters for the Ben R euven line shape come from 
IDeboer an d Stcffcs (|1994j ). We choose as our nominal 
atmospheric model no H2S absorption, and additionally 
investigate models with absorption due to 10 and 50 
times solar H2S abundances. We note that high H2S 
abundances give a best fit to the cm wavelength region 
(Section EU Fig. HJ. 

3.3. Thermal profile 

The observed shapes of the CO absorption lines de- 
pend strongly on both the CO abundance and the 
temperature-pressure (TP) profile in the atmosphere. 
Therefore, we look at a range of TP profiles to under- 
stand the effect of uncertainties in the TP profile on our 
derived CO altitude profile. Our nominal TP profile is 
th at of IFletcher et all (|2010f ): they adopt the TP profile 
of iMoses et al.l (|2005l ) below the tropopause and retrieve 
the stratospheric temperature profile from mid-infrared 
AKARI spectra. Fits to HD lines in Neptune's 51-220 
/zm spectrum from Herchel/PACS (|Lellouch et alJl2010( ) 
are sensitive to the 10-500 mbar levels; these authors 
favor a thermal profile much like that of IFletcher et al.l 
(|2010l) in the stratosphere as well. 

Figure [5] shows a selection of published TP profiles, in- 
cluding our nominal profile and the profi l es used in the 
recent CO studies oflLellouch et al.l (|2005| ) , iMarten et al.1 
(|2005f l and lHesman et all (|2007f ). Published profil es gen- 
erally match the Voyager RSS occultation profile (|Lindall 
1992) at ~1 bar (T = 71.5 K), and assume an adia- 
batic extrapolation down to deeper leve ls. Ou r nom - 
inal profile is slightly warmer than the iLindall (|1992f ) 
profile at the tropopause, which could be due to sea- 
sonal change s in the atmospheric temperature since the 
Voya ger era ()Orton et al.ll2007at lHammel and Lockwoodl 
2007). At altitudes above 1 bar, the published profiles 
diverge greatly. We consider four additional TP profiles 
that span the range of published profiles; each of these 
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Fig. 4. — Model atmosphere calculations. All models include CO and H2 collision-induced absorption (see Section [3]2}; the red 
line represents these two opacity sources only. Opacity due to a 10 and 50 times protosolar abundance of H2S (cyan, dashed 
and solid lines, respectively); a 50 times protosolar abundance of H2O (orange); a protosolar abundance of NH3 (green) and 
a 0.1 times protosolar abundance of PH3 (blue) are included. The black dashed model contains all of these absorptions (note 
that a 50 times protosolar abundance was used for H2S, rather than the 10 times sola r value). A compilatio n of microwave 
data taken before 1991 is plotted for comparison (diamonds); these data are described in lde Pater et all (|1991aD and references 
therein. The right plot zooms in on millimeter wavelengths, and shows the importance of H2S opacity at these wavelengths. 



test profiles matches the nominal profile at altitudes be- 
low 1 bar. The 'extreme low' and 'extreme high' profiles 
decrease/increase steadily to ±25 K from the nominal 
profile, respectively, between 1 and 0.003 bar and then 
remain at ±25 K from nominal as the pressure continues 
to decrease. The 'moderate low' and 'moderate high' pro- 
files decrease/increase steadily to ±10K from the nom- 
inal profile, respectively, between 1 and 0.008 bar and 
then remain at ±10 K from nominal as the pressure con- 
tinues to decrease. These offsets are illustrated in Fig. 
El While these test profiles envelope the published pro- 
files and their error bars, they do not represent all of the 
possible TP p rofile shapes. For example, the profile of 
IBezard et all (|1999ft is ~20 K cooler than our nominal 
profile around 10 -3 bar, and warmer than our nominal 
profile at higher and lower pressures. The provisional 
temperature profile from the analysis of 2005 Spi t zer In - 
frared Spectrom eter data, rep o rted i n lLine et al.l ((2008), 
is similar to the IBezard et al.l (| 1999T ) profile. We do not 
test all of the published TP profiles, but we do consider 
the IBezard et alJ {T999) profile in addition to our nom- 
inal profile an d four test profiles. T he p rofiles used in 
the an alyses of iLellouch et al.l (|2005| ) and iHesman et all 
(|2007f ) are most similar to our 'moderate cool' profile. 

At pressures greater than 1 bar, where the temperature 
profile is unknown, we assume the atmosphere is con- 
vective and extrapolate adiabatically. We consider two 
poss ibilities: a dry a diabat (nom inal) and a wet adiabat 
(see Ide Pater et ail (fl991a. 2005) for a detailed descrip- 



tion). In the adiabatic extrapolation, we assume that the 
specific heat of hydrogen is near that of normal hydrogen 
(though the opacity for equilibrium h ydrogen is used- this 
is the "intermediate" hydrogen ca se (|Massie and Huntenl 
fl982T : lde Pater and MitchellT^l ). 

3.4. Disk averaging 

To account for variations in viewing geometry on Nep- 
tune's disk-averaged spectrum, we calculate Tb{v,[i) for 
25 different viewing angles /i, which represent the aver- 
age spectra within 25 concentric rings. Doppler broad- 
ening due to the planet's rotation, which has a velocity 
of V eq = 2.7 km s _1 at Neptune's equator, affects the 
shape of the e mission peak in the ce nter 2 MHz of the 
line. Following iMoreno et a l. (2001), the velocity at a 
given distance x from the central meridian is: 

V = ^V eq (3) 

where 

x = Rcos(lat) sm(Along) cos(SEPL) 

We define R as the planetary radius, lot as planetary lat- 
itude, Along as longitude from the central meridian and 
SEPL as the latitude of the sub-earth point. We divide 
the disk into 100 values of x. At each x, we calculate the 
average spectrum given the relative contribution of each 
representative viewing angle fj,. We then shift the spec- 
trum in frequency given the value of V(x). Finally we 
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Fig. 5. — Selected temperature- pressure (TP) profi les for 
Neptune from the literature. The iMarten et all (|2005f ) ther- 
mal profile (teal) was derived from th eir submillimet e r data 
concurrentl y with CO abunda nce. ILellouch et al.l (2005) 
adopted the lBezard et al.l ([19911 ) T P profile (orange) in their 
study; the the rmal profile use d bv [Hesman "e t al. (2003) was 
based on the iBurgdorf et alJ l|2003l) profile (brown). Our 
nominal choice for this work is the lFletcher et all (120101) pro- 
file (t hick black line), and we also consider the Be zard et alj 
(1999) profile (thick red line). At pressures greater than 1 
bar, we follow a dry adiabat (black, short-dashed line). 

coadd the Doppler-shifted spectra to get our final disk- 
averaged spectrum. Once the disk-averaged spectrum is 
computed, models are convolved to the instrumental res- 
olution (^1 MHz near the peak, ^33 MHz in the wings). 

4. ANALYSIS 

To determine Neptune's vertical CO profile, we per- 
form robust non-linear least-squares fitting to the data 
weigh ted by the data er rors using the MPFIT IDL pack- 
age (jMarkwardtl [2009). We first consider models in 
which the CO mole fraction tico/^tot * s ne ld constant 
throughout the atmosphere. This corresponds to the case 
where Neptune's observed CO primarily comes from ver- 
tical mixing from deep levels in the atmosphere. In ad- 
dition to the value of nco/^TOT, these fits have two ad- 
ditional free parameters. These are amplitude correction 
factors for the 1- and 3-mm data, which account for the 
uncertainty in the absolute flux calibration of the data. 
In practice, we find the values of these flux density cor- 
rection parameters to be in the range of 5-10% at both 
wavelengths. 

The least-squares fitting procedure is performed for 
models that use our nominal TP profile as well as each 
of the test thermal profiles described in Section 13.31 We 
also perform the fit on the 1-mm and 3-mm datasets 
separately. We report the fitted values for the CO mole 
fraction, 1- and 3-mm amplitude correction factors, and 
the statistical errors from the covariance matrices of the 
fits in Table [3] We also report the reduced chi-squared 



Fig. 6. — Difference between selected TP profiles and our 
nominal TP profile. Line colors are the same as in Fig. [5] 
several additional thermal profiles are also shown. Published 
alternate profiles/uncertainties are included as broken lines. 
The light grey lines indicate our nominal test profiles, chosen 
to span the range of published profiles above 1 bar (see Section 
[3]for details). From left to right, these are the 'extreme cool', 
'moderate cool', 'moderate warm' and 'extreme warm' TP 
profiles. 

(x 2 ) of the best fit: 

X 2 = — V Vm (4) 

* M-N ^ a 2 V ' 

m _Q meas,m 

where M is the x 2 number of data points, and N is the 
number of fit parameters (in this case, three), so that 
M — N is the number of degrees of freedom (DOF) for 
the fit. The parameter 6y m is the difference between the 
data and model value for point m, and a meas .m is the 

measurement error for point m. Our values of x 2 are 
greater than one: this is discussed in detail in Sec t ion [51 
Following previo us authors ([Lellouch et al.l 120051 : 
iHesman et al.l I2007T ) we repeat the fitting process using 
a two-level CO profile. In this case the CO abundance is 
assumed to be constant within each of two levels, above 
and below some transition pressure. Fits of this form 
are representative of the scenario in which CO is pro- 
duced in the stratosphere and diffuses downward, along 
with being mixed upward from the interior. In addition 
to the CO mole fractions and the 1- and 3- mm am- 
plitude correction factors, the pressure that defines the 
transition between the two atmospheric levels must also 
be determined. We find that our fitting procedure does 
not perform well when we allow this transition pressure 
to be a free parameter; the fitting program will rarely 
vary the pressure level from its start value. Therefore, 
for each model thermal profile, we perform a series of 
fits: in each run we fix the pressure level, testing transi- 
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TABLE 3 

Best-fit one- level CO profiles. The CO mole fraction is held constant as a function of altitude (Section 16. ip . The 1- and 3-mm 
flux density factors are the factors by which the data must be multiplied in the best-fit solution. These factors account for 

errors in the gain calibration of the data, which may be as much as 20% (see Section [2J. x 2 1S the reduced \ 2 °f the best-fit 
model, as defined in Eq. (4). The thermal profile used in each model is specified. 



Data 


TP profile 


CO (ppm) 


1-mm flux density 
factor 


3-mm flux 
factor 


density \ 2 


1 mm 














nominal 


0.50 ±0.02 


1.050 ±0.003 




15.3 


3 mm 














nominal 


0.37 ±0.05 




1.09 ±0.004 


14.6 


1 mm - 


- 3 mm 












nominal 


0.49 ±0.02 


1.052 ± 0.003 


1.081 ± 0.002 


15.3 




Bezard et al. f 19991 


0.45 ±0.01 


1.055 ±0.003 


1.083 ±0.002 


15.4 




extreme cool 


0.264 ± 0.008 


1.061 ±0.003 


1.090 ± 0.002 


18.3 




moderate cool 


0.37 ±0.01 


1.055 ±0.003 


1.085 ± 0.002 


16.2 




moderate warm 


0.60 ± 0.02 


1.056 ±0.003 


1.079 ± 0.002 


16.5 




extreme warm 


0.67 ±0.03 


1.069 ±0.004 


1.080 ± 0.003 


21.5 



tion pressures in steps of 0.1 in logP(bar) from -3.0 (1 
mbar) to 0.0 (1 bar). The best-fit solutions as a function 
of transition pressure are shown in Fig. [7] for a selec- 
tion of thermal profiles. We find that for a given TP 

profile, the x 2 value has a global minimum, and all the 
parameters are smooth functions of the log P(bar) of the 
transition pressure. We therefore take the solution with 

the lowest x 2 value as our overall best-fit model for each 
test TP profile; for example, for our nominal profile (no 

H2S), the minimum x 2 reveals a transition pressure of 
logP(bar) = 0.9, or 0.13 bar (Fig. top panel). This 
corresponds to a CO abundance of 1.1 ppm at altitudes 
above 0.13 bar (middle panel) and 0.1 ppm at altitudes 
below 0.13 bar (lower panel). These best-fit values are 
reported in Table |4j As when modeling a constant CO 
mole fraction, we repeat the two-level fit for the 1- and 
3-mm datasets separately. We also test the effect of H 2 S 
absorption on the best-fit two-level CO profile, by repeat- 
ing the two-level fits assuming a 10 and 50 times solar 
abundance of H2S. These results are reported in Table 0] 
as well. 

Finally, we use a model to derive several 'physical' CO 
profiles. To model the vertical diffu sion ra t e, the eddy 
diffusion coefficient {K ) profiles of IMosesI (|1992l ) and 
IRomani et al.l ()1993D . which were developed to match 
photochemical models of observed hydrocarbon distribu- 
tions, ajejis^d_(Fi_g_ 1 j8|^Mokcular diffusion is included 
using iMarrero and Masonl (|1972D : 




-2.5 -2.0 -1.5 -1.0 

Transition pressure (log P(bar)) 



D C o = 
1 



1.539 x io- 2 T 1548 



P{atm) (ln[T/3.16 x 10 7 ]) 2 e -2.8/T e i067/T^ 



cm 2 s 1 



(5) 



Using these profiles, we solve the diffusion equation: 



, T r ( dnco 



D 



CO 



-"-co 

dn C Q 
dz 



1 

H 



nco 



H, 



1 dT 
f~dz_ 
1 1 dT 
^ T~dz~ 



CO 



(6) 



where nco is the number density of CO molecules in 
cm~ 3 , z is the altitude in cm, H is the total atmo- 



Fig. 7. — Best-fit parameters as a function of transition pres- 
sure level for the nominal thermal profile (black), both with 
(dashed) and without (solid) H2S opacity; and for the extreme 
warm (red) and extreme cool (blue) thermal profiles. Pres- 
sure is given in units of log(P(bar)), from log(P(bar)) = —3.0 
(1 mbar) to log(P(bar)) = 0.0 (1 bar). 

spheric pressure scale height, and Pco is the scale height 
of CO. We then find the best-fit value for the exter- 
nal flux 4>co from the upper atmosphere, allowing for 
an additional CO contribution from the planet's interior 
( 71 Co/ 77 -tot|( z =o))- Sample physical CO profiles, assum- 
ing an influx rate of <f>co = 10 8 cm~ 2 s -1 and no internal 
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TABLE 4 

Best-fit two-level CO profiles. The CO mole fraction transitions between a high altitude value ("CO above") and a deep 
atmosphere value ("CO below") at the best-fit transition pressure P which is given in units of logP(bar) as well as in bars 
(Section 16.216.31) . As in Table [3] the 1- and 3-mm flux density factors are the factors by which the data must be multiplied to 

match the model in the best-fit solution, and \ 2 ls the reduced \ 2 of the best-fit model, as defined in Eq. (4). 



Data 




TP profile 


CO below (ppm) 


CO above (ppm) 


log P(bar) 


P(bar) 


1-mm flux density 
factor 


3-mm flux density 
factor 


X 2 


1 mm 
























nominal 


0.09 ± 0.03 


1.12 ± .04 


-0.9 


0.13 


1.079 ± 0.003 




10.7 


3 mm 
























nominal 


0.00 ± 0.01 


0.89 ± .08 


-1.0 


0.100 




1.100 ± .0002 


9.3 


all 




nominal 


0.08 ± 0.03 


1.06 ± 0.04 


-0.9 


0.13 


1.082 ± 0.003 


1.088 ± 0.002 


11.1 






Bezard et al. ( 1999) 


0.12 ± 0.03 


1.05 ± 0.05 


-1.0 


0.10 


1.080 ± 0.003 


1.088 ± 0.002 


11.4 






extreme cool 


0.285 ± 0.007 


1.9 ± 0.1 


-2.6 


0.0025 


1.062 ± 0.003 


1.089 ± 0.002 


12.3 






moderate cool 


0.30 ± 0.01 


1.37 ± 0.09 


-1.6 


0.025 


1.066 ± 0.003 


1.085 ± 0.002 


12.7 






moderate warm 


0.001 ± 0.03 


1.13 ± 0.04 


-0.7 


0.20 


1.094 ± 0.003 


1.090 ± 0.002 


12.3 






extreme warm 


0.00 ± 


0.93 ± 0.03 


-0.4 


0.40 


1.100 ± 0.003 


1.092 ± 0.002 


18.7 


+ 10x solar 
+ 50 X solar 


H 2 S a 
H 2 sE 


nominal 
nominal 


0.00 ± 
0.00 ± 


1.08 ± 0.02 
1.08 ± 0.02 


-0.8 
-0.8 


0.16 
0.16 


1.036 ± 0.002 
1.039 ± 0.002 


0.904 ± 0.001 
0.909 ± 0.001 


9.3 
9.3 



a Opacity due to 10 X solar enrichment in H2S included 
b Opacity due to 50x solar enrichment in H2S included 
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Fig . 8. — Eddy diffusion coe fficient profiles from IMosed 
(1992]) and lRomani et alj (fl993h . Also included is one of our 
'test profiles', designed to produce a CO profile more similar 
to our best-fit two-level CO profile (see Section 16.41 for dis- 
cussion). The tropopause pressure values of 0.11/0.20 bar are 
indicated. 

contribution of CO, are illustrated in Fig. [9] To demon- 
strate the effect of internal CO on these physical CO pro- 
files, we show the same profiles with 0.2 ppm of internal 
CO added for two of the eddy profile cases. Best-fit val- 
ues of <f>co and nco/«TOT|( z =o) from our least-squares 
modeling are summarized in Table [5] 

5. ERRORS AND UNCERTAINTY 

As described in Section [2j the least-squares fitting uti- 
lizes the residual map rms scatter to weight the data for 

calculating \ 2 . This weighting scheme compensates for 
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Fig. 9. — Sample CO profiles that result from eddy profiles 
in Fig. [8] assuming <f>co = 3 x 10 8 cm"^" 1 and no internal 
contribution to the CO. To illustrate the effect of internal 
CO on the CO profile, we also show the CO profile for the 
Romani B eddy profile (solid red line) with the addition of 
0.2 ppm of internal CO (solid grey line); and the 'test eddy 
profile' (orange dot-dashed line) with 0.2 ppm of internal CO 
included (grey dot-dashed line). The best- fit parameters for 
physical CO profiles depend on the TP profile in addition to 
the eddy profile; fits are summarized in Table[5] For reference, 
the best-fit two-level profiles from our nominal (no HhS, cyan 
line) and extreme cool (blue line) TP profile fits are shown. 

the difference in signal-to-noise between the narrow (~ 1 
MHz) and wide (~ 33 MHz) frequency channels, which 
we expect to differ by a factor of ^6. Noisy channels 
that are not flagged during data reduction are down- 
weighted. Day-to-day variations due to weather or array 



12 



S. H. Luszcz-Cook and I. de Pater 



TABLE 5 

Best-fit physical CO profiles. The external flux </>co and the mole fraction of CO brought up from the deep atmosphere 
("internal CO") that produce the best match between the spectrum and model are given, for a selection of thermal and eddy 

diffusion coefficient profiles (see Section PS . 4 p . \ 2 is the reduced x 2 of the best-fit model, as defined in Eq. (4). The flux density 

scale factors for these fits are not shown. 
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performance, and differences in the noise between 1 and 
3 mm, are also accounted for. However, we expect that 
these weights underestimate the total data uncertainties. 
While the overall calibration offset is a free parameter in 
the least-squares fits, additional systematic errors in the 
flux densities, for example due to errors in the band- 
pass calibration, could affect the relative flux between 
frequency windows. A visual comparison of individual 
1-mm tracks (see Fig. [3] ) shows small day-to-day offsets 
in the average flux density in a given frequency window, 
for example at 230 GHz and 235 GHz, that are likely due 
to this type of systematic error. 

Additionally, gain calibration offsets between each of 
our 1-mm tracks could also cause systematic, frequency- 
dependent errors in the final spectrum. Since we expect 
a residual uncertainty of 3% in the relative gain calibra- 
tion between our 1-mm datasets, we tested the effect of 
such errors using a Monte Carlo method. We assume 3% 
random errors in the relative gains between data sets and 
perform 300 trials of our constant- value nco/^TOT fit. 
We find that the fit solution is stable to the inclusion of 
these errors, but the statistical uncertainty of the best- 
fit CO abundance increases by a factor of 2. While this 
experiment is too time intensive to repeat for the other 
model cases, this result suggests that 1) in general we 
are underestimating the errors in the fits by ignoring the 
error in the relative gains between 1-mm tracks, and 2) 
at least in this case, including the errors in the relative 
gains between 1-mm tracks leads to an increase in the 
uncertainties, but does not change the best-fit solution. 

Because of the systematic errors, we consider our val- 
ues of <Jmeas,m to be weights rather than true estimates of 
the uncertainties in our data points. This is reflected in 

the values of x 2 for our fits, which we would expect to be 
1.0 for a model that fully captures the data if the values 
of tJmeas.m were reflective of the true data uncertainties. 

A x 2 of 11.1 for our best-fit, two- level model implies that 
if this model were a perfect representation of the data, 



then the data weights (J m eas,m underestimate the total 
data uncertainty by a factor of 3.3. However, the high 

value of x 2 1S likely due to deficiencies in the model in 
addition to an underestimate in the data errors. Factors 
that contribute to this include errors in the continuum 
opacity: either in the mixing ratios of the continuum 
species or in the coefficients/line profiles used to derive 
continuum absorption. Our tests with H2S, a species 
which we observe to affect our millimeter spectrum but 
whose abundance is poorly constrained, provide one ex- 
ample of this. Errors in the shape of the TP profile and 
CO profile can also contribute, since we only investigate 
a finite number of TP and CO profile shapes. 

In summary, the x 2 parameter does not reach the the- 
oretical expectation of x 2 ~ 1 for a perfect fit because 1) 
while the values of cr meas m adopted from the rms scatter 
in the maps are good estimates of the relative dispersions 
of the data points, they do not characterize the total in- 
trinsic uncertainties; and 2) the models (even our best 

models) are not perfect. Regardless, x 2 is a useful mea- 
sure of the variance of the data point residuals (Sy^, see 
above) for our different fits. We encourage the reader 

to use both the reported x 2 values and the plots when 
evaluating the success of various models in reproducing 
the data. 

6. RESULTS 

The results of our model fits, which are summarized in 
Tables 3-6, are illustrated in Figs. 10-21 and discussed 
below. Plots of the CO spectra show both the unbinned 
single-channel data with uncertainties cr meaStm , and the 
data binned to 3 times the channel width; 1- and 3-mm 
flux densities are typically scaled by the amplitude cor- 
rection factors from the nominal TP profile best two-level 
fit (not including H 2 S), before being converted to bright- 
ness temperature. For ease of comparison, additional 
model spectra are also scaled to align with this nominal 
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Fig. 10. — Models of the millimeter spectrum corresponding to the best-fit constant (0.5 ppm) CO profile in blue; and two-level 
(1.1 ppm above the 0.13 bar level, 0.1 ppm deeper than 0.13 bar) CO profile in red. EhS absorption is omitted here. For plotting 
purposes, the 1- and 3-mm data were scaled up by 8.2% and 8.8% respectively, to match nominal TP, best-fit two-level model; 
likewise, subsequent models are all scaled by the relative flux correction factors to permit comparison. This plot shows that a 
two-level profile fits the data significantly better than a one-level profile. 



two-level model. 

6.1. Constant CO models 

Our best-fit model for a constant CO distribution uses 
the nominal TP profile and a CO mole fraction of 0.5 
ppm. The best-fit CO abundance decreases for cooler 
TP profiles, and increases for warmer TP profiles, within 
a range of 0.3 — 0.7 ppm (Table |3]). Figure fTOl illustrates 
that a constant CO profile that has the appropriate line 
depth is generally too broad in the wings. The emission 
peak at line center is also a poor match to the data. 

6.2. Two-level CO models, no H2S 

Model spectra produced with our two-level CO pro- 
files offer an improved fit to the data; this is also il- 
lustrated in Figure [10] Overall, our nominal TP pro- 
file allows the best fit to the data, in a v 2 sense, with 



the iBezard et al.l (|1999t ) TP profile solution a close sec- 
ond. The best two- level model solution has a x 2 of 11.1, 
with 1.1 ppm of CO in the upper atmosphere and 0.1 
ppm CO in the lower atmosphere. The transition pres- 
sure for this model is located near the tropopause, at 
\ogP(bar) = -0.9 (P = 0.13 bar). The best-fit parame- 
ters for several transition pressures near the x 2 minimum 
value of log P(bar) = —0.9 (P — 0.13 bar) are presented 
in Table [5] the model CO (2-1) spectrum for three of 
these solutions is plotted against the scaled l-mm data 
in Fig. [HI The good match of all of these spectra to the 
data implies that there is a range of transition pressure 
levels that produce quality fits to the data. This is also 

clear from the best-fit x 2 values; for example, the fit so- 
lution at log P(bar) — —0.8 (P = 0.16 bar) also has a x 2 
of 11.1, with a slightly lower best-fit CO abundance of 
1.0 ppm in the upper atmosphere and 0.0 ppm CO in the 
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TABLE 6 

Best-fit two-level CO profiles for a series of transition pressure levels, for models using our nominal thermal profile 
(|Fletcher et al. 2010). The CO mole fraction transitions between a high-altitude mole fraction ("CO above") and a 
deep-atmosphere mole fraction ("CO below") at a transition pressure P which is given in units of logP(bar) as well as in bars 
(Section 16. 2[) . Transition pressures in the range of logP(bar) = —3.0 to 0.0 were tested with the model, and a subset of these 
tests is shown here. As in Table [3j the 1- and 3-mm flux density factors are the factors by which the data must be multiplied 

to match the model in the best-fit solution, and x 2 ls t ne reduced \ 2 °f l ne best-fit model, as defined in Eq. (4). 
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0.27 ±0.02 


1.27 ±0.07 


1.074 ± 0.003 
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lower atmosphere. Therefore, the best-fit answer alone 
does not fully characterize the range of solutions allowed 
by our data. We also find that models using our nominal 

thermal profile produce solutions with lower values of x 2 
than the best solutions from models using other thermal 
profiles, for a range of transition pressures (see Tables [4] 
and[6]). These solutions have CO mole fractions of 0.0-0.3 
ppm in the troposphere and 0.8-1.3 ppm in the strato- 
sphere, with level transitions between 0.06 and 0.25 bar. 

This very shallow minimum in x 2 as a function of transi- 
tion pressure is likely a result of the thermal profile being 
nearly isothermal in this region around the tropopause. 
In addition to the default assumption of a dry adiabat, 
we also test thermal profiles that used a wet adiabat to 
extrapolate to high pressures. We find that, in general, 
using a wet adiabat for the thermal profile does not alter 

the models significantly, and increases the values of x 2 
slightly. 

Of all of our test thermal profiles, models using the 
extreme warm profile agree least well with the data, in a 

X 2 sense. This thermal profile is also the least consistent 
with the published TP profiles (see Fig. [6]). The best-fit 
transition pressure for this thermal profile is at 0.4 bar, 
which is below the tropopause, and therefore probably 
unphysical. The best-fit moderate warm CO profile is 
qualitatively similar to the best model from the nominal 
TP profile, with 1.1 ppm of CO in the stratosphere and 
0.0 ppm CO below the tropopause. Figures Q2] and [13] 
show that both warm TP profiles produce line profiles 
that are too high in the near wings of the line (< 200 
MHz from center), particularly at 1 mm. The extreme 
warm profile is also too broad in the far wings of the CO 
line. 

Model fits using the cooler TP profiles favor two-level 
profiles with level transitions occurring well above the 
tropopause; the transition for the best-fit, moderate cool 
thermal profile case is at 25 mbar; and at 2.5 mbar for 
the extreme cool profile. Fits using cooler TP profiles 
favor higher CO mole fractions in both levels of the at- 
mosphere: the best-fit solutions give 0.3 ppm CO in the 
lower level of the atmosphere for both thermal profiles, 
and 1.4 and 1.9 ppm CO in the upper atmosphere for 
the moderate and extreme TP profiles, respectively. As 
Figs. [12] and [13] show, models using the cooler TP pro- 
files have the correct absorption depth and match the 
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Fig. 11. — Spectral models for transition pressures of 
log(P(bar)) = -0.8, -0.9 and -1.0 (0.16, 0.13, and 0.10 bar, 
respectively), and our nominal TP profile, shown for the CO 
(2-1) line (see Table [6]). Solutions for log P(bar) = -0.8 
(green dashed line) and —0.9 (solid magenta line) are some- 
what more favorable in the far wings (2-3 GHz from line cen- 
ter), but all three models match the data quite well. 



emission peak well; but are too broad in the far wings, 
particularly at 1 mm. 

Given these findin gs, we conclude that the data are 
best matched by the lFletcher et all (]2010D thermal pro- 
file, with O.l^o; 2 ppm of CO in the troposphere and 
l.lio'3 PP m °f CO in the stratosphere, as indicated by 
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Fig. 12. — Comparison of best-fit two-level models to the CO 
(2-1) data using different TP profiles. 

the range of values in Table [6j Fits to the 1-mm and 3- 
mm lines separately, which are performed using the nom- 
inal TP profile, give consistent best-fit solutions. Visual 
inspection of all of the models generally indicates bet- 
ter agreement with the data at 1 mm; however since the 

3-mm data have lower signal-to-noise, the \ 2 value for 
the 3-mm data alone is actually lower than for the 1-mm 
data. Best-fit corrections to the flux density are 6-10% 
at 1 mm and 9-10% at 3 mm (Table gj. 

6.3. Two-level models, H2S included 

As discussed in Section|31 we test the effect of including 
opacity due to a 10 and 50 times protosolar enrichment 
in H2S. We find that the main effect of adding H2S on 
the best-fit model parameters is to change the amplitude 
scaling factors for the 1- and 3-mm data. Most notably, 
without H2S, the 3-mm data are lower than the model 
by ~ 9%, whereas with 10 x protosolar H2S the data 
are higher than the model by ~ 9% (Fig. H"4"|) . Abun- 
dances of H2S higher than 10 x the protosolar value do 
not further change the millimeter spectrum, due to con- 
dens ation into the H 2 S ice cloud at pressures of a few 
bars (|de Pater et al.lll991alh 

Despite the significant effect of H2S absorption on the 
3-mm continuum, we find that the best-fit CO solution is 
not strongly affected by the inclusion of H2S absorption: 
the best-fit CO vertical profile has 1.1 ppm CO in the 
stratosphere and 0.0 ppm of CO in the troposphere (Ta- 
ble SJ. The best- fit CO model spectra, scaled to match 



Fig. 13. — Comparison of best-fit two-level models to the CO 
(1-0) data using different TP profiles. 
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Fig. 14. — Best-fit solutions for models containing no H2S 
absorption (solid cyan), 10 times solar H2S (blue dashed) and 
50 times solar H2S (red dotted). In this Figure, the data and 
models have not been scaled by any flux density correction 
factor. Generally, all models are higher than our raw data at 
1 mm. At 3 mm, the no-HbS model is higher than the raw 
data, whereas the HbS-enriched models are lower. The 10 and 
50 times protosolar H2S models match one another closely. 

the best two- level fit for the no-H 2 S case, are plotted 

with the (scaled) data in Figs. [TBI and UM The x 2 value 
for the fit improves with the addition of H2S, which is 
primarily due to an improved fit to the shape of the far 
wings of the CO (1-0) line (Fig. ITS"]) . There is also some 
improvement to the fit at the center of the CO (1-0) line. 

6.4. Physical models 
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Fig. 15.— Model fits to the CO (2-1) line as in Fig. HH 
except that the data and HbS-enriched models are scaled 
to match the brightness temperature of the best-fit, no-IUS 
model. 

Physical C O pr ofiles based on the d iffusion models of 
iMosed (|1992D and lRomani et all (|1993f ) produce the best 
model spectra when cool thermal profiles are used. The 
reason for this can be inferred from Fig. |9j for all of these 
diffusion models, the mole fraction of CO in the atmo- 
sphere from an external source falls off at pressures less 
than ^10 mbar, which is much higher in altitude than the 
transition level for the best-fit two-level model using our 
nominal thermal profile, but is consistent with the transi- 
tion levels for the best-fit moderate and extreme cool TP 
profile models. The internal CO contribution for these 
fits is roughly independent of the eddy diffusion profile, 
with mole fractions of 0.3 and 0.4 ppm for the extreme 
and moderate cool thermal profiles, respectively. These 
values are very similar to the CO abundances found for 
the two-level fits using the same thermal profiles. The ex- 
ternal CO flux is dependent on the choices for the eddy 
diffusion coefficient and thermal profiles: values range 
from 0co = 0.5 - 7 x 10 8 cm" 2 s _1 . 

To produce a physical CO profile that is more like our 
best-fit two-level CO distribution, we created two 'test' 
eddy diffusion coefficient profiles (Fig. [5J. The pro- 
files are designed to have a more rapid increase in the 
mixing from the diffusion rate minimum just above the 
tropopause, to allow high concentrations of externally 
supplied CO to reach deeper levels. No attempt is made 
to evaluate the physical likelihood of such an eddy profile. 
The two test profiles are identical except for the location 



Fig. 16.— Same as Fig. [151 except for the CO (1-0) line. 

of the tropopause, which is defined by a sharp transition 
from fast to slow mixing. In addition to using a value 
for the tropopause of 0. 11 bar as in the IMosesI Jl992) 
and iRomani et ail (|1993j ) cases, we also try locating the 
tropopause level at 0.20 bar. Using these test eddy pro- 
files, we are able to produce physical CO models using 
the nominal TP profile that are in better agreement with 
the data. The best-fit parameters for the test profiles are 
<j> = 10 - 20 x 10 s cm^V 1 and CO mt = 0.2 - 0.3 ppm. 
Unsurprisingly, the test eddy profiles, which are designed 
to be used with the nominal TP profile, do not do as well 
when used with the cool TP profiles. A selection of the 
best model spectra based on physical CO distributions is 
shown in Figs. [T7l and fT8l 

7. DISCUSSION 

Our analysis of the CARM A CO (2-1) and (1-0 ) data 
indicates a preference for the iFletcher et al.l (|2010t ) ther- 
mal profile in Neptune's upper atmosphere over warmer 
and cooler profiles. Good fits are characterized by CO 
mole fractions of 0.0-0.3 ppm below the tropopause and 
0.8-1.3 ppm in the stratosphere. This is true for the in- 
dependent fits of t he 1- and 3-mm line s, as well as the 
combined fit. The Bezard et al.l (|1999[ ) thermal profile 
and our moderate warm profile give best-fit CO profiles 
that are consistent with the nominal thermal profile re- 
sults. Cooler thermal profiles also produce acceptable 
fits to the data, but models using the cooler TP profiles 
are generally too broad in the far wings of the CO lines. 
The vertical CO profiles derived using our moderate and 
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Fig. 17. — Model fits to the CO (2-1) line from a selection 
of best-fit physical CO profiles (see Table [5]). 

extreme cool thermal profiles are qualitatively different 
than the nominal TP profile best-fit solution: we find 
CO mole fractions of 0.3-0.4 ppm in the troposphere and 
lower stratosphere, and 1.4-1.8 ppm above 2.5-25 mbar. 
Using a physical diffusion model, we determine that the 
stratospheric CO abundances for all thermal profiles cor- 
respond to an external CO flux of c/>co = 0.5 — 20 x 10 8 
cm _2 s . Internal CO contributions for the physical CO 
profile solutions are typically 0.2-0.4 ppm. 

Plots of the 3-mm spectrum (e.g. Figs. [TUl Q2] [TBI 
[T8|) show greater deviations between the data and models 
than we find at 1 mm. We attribute this to the fact that, 
at 3 mm, we probe deeper levels of the atmosphere (Fig. 
[JJ), where the atmospheric opacity is not well constrained. 
However, we are encouraged by the fact that the indepen- 
dent fits to our two CO lines give consistent CO profile 
solutions. Additionally, we find that the inclusion of 10- 

50 x solar H 2 S opacity improves the x 2 goodness-of-fit, 
but does not significantly affect the best-fit CO vertical 
profile. 

7.1. Comparison with previous results 

In Fig. [T9]we compare our best-fit two-level CO profiles 
with previously reported values and profiles. Figures [50] 
and [5T] use several of these recent results for the CO 
profile in conjunction with our radiative transfer code to 
produce model spectra. In each case we fit only the data 

amplitude scale factors at 1 and 3 mm; the lowest \ 2 f° r 
these literature CO profiles and our data are presented 



Fig. 18.— Same as Fig. [T71 except for the CO (1-0) line, 
in Table El 

Of the three most recent pre vious measur e ments 
of millimeter CO li ne shapes dMarten et al.l 120051 : 
Lellouch et all 120051: iHesman et al.l 120071 ). only 
Mar ten et al. I (|2005F found a result consistent with 
a constant CO mole fraction. However, their data only 
cover wavelengths up to 50 MHz away from the center 
of the CO (4-3) line. From the CO line contribution 
functions of the first three rotational transitions (Fig. 
UJ), we estimate that their data are only sensitive down 
to pressures of a few tens of mbar, which implies that 
they would not detect a troposphcric decrease in the 
CO abundance. Our best-fit two-level, nominal TP 
profile model has an abundance of 1.1 pp m in the upper 
atmos phere, which agrees well with the iMarten et al.l 
(2005) value. As Figs. HOI and CUJ and Table illustrate, 
this profile is the least successful of those tested at 
reproducing our data. 

Our two-level solutions, which have a smaller abun- 
dance of CO in the lower atmosphere, are qualitatively 
consistent wit h the results o f the two r ecent millime- 
ter stu dies by iLellouch et al.l (|2005t ) and IHesman et al.l 
(|2007| ). Quantitatively, there are important differences 
between t hese two previo u s wor ks and with our own re- 
sult. The IHesman et al.l (|2007D best-fit solution has a 
str atospheric CO abund ance that is inconsistent with 
the ILellouch et all (|2005t ) result, as well as a lower pres- 
sure that defines the transition between the upper and 
lower levels of a two-level model. Our analysis demon- 
strates that there is a strong relationship between the 
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Fig. 19. — Selected published CO profiles with errors (left) , compared to our b est-fit profiles (right). Grey boxes indicate 
all the pu blished results from prior to 1996, as descri bed in [C ourtin ct al. (1996). Colored regio ns indicate the best-f it CO 
profiles of Mar ten et al.l (2005) (green, cross-hatched) , [Lellouch et al.l (|2005) (orange, shade d) and IHesman et all (|2007D (ma- 
genta, hatched ) ; alon g with their uncertainties. The symbols indicate the values reported bv lFletcher et al.l (|2010f ) (stars) and 
ILellouch et al.l (2010) (diamond) for the stratospheric CO abundance. On the right, our best-fit two-level profiles are shown for 
several different thermal profiles, as well as for the case where absorption due to 50 x solar H2S is included. 



stratospheric CO mole fraction and the pressure that de- 
fines the transition between the upper and lower levels 
of a two-level model. Furthermore, the goodness-of-fit, 

as measured by x 2 , tends to be relatively constant over 
a wide range of transition pressures, so that small varia- 
tions in the data may lead to significant changes in the 
derived CO abundances. This leads to uncertainties in 
the CO abundances that are greater than the 20% uncer- 
tainties adopted by ILellouch et all (|2005). and may be a 
do minant source o f inco nsiste ncy between the sol utions 



of ILellouch et al.l (J2005) and IHesman et al 



comparison w i th ou r results, ILellouch et al 



(2007). In 
2005) and 



IHesman et al.l ([2007D have higher tropospheric CO abun- 
dances (0.5 ±0.1 and 0.6 ± 0.4 ppm, respectively) and 
transition altitudes (20 and 6 mbar, respectively) than 
our nominal best- fit solution of less than 0.1 ppm in the 
troposphere, with a transition near the tropopause. In- 
terestingly, our best-fit solution using a moderate cool 
th ermal profile, which is simi lar to the TP profile s used 
bv ILellouch et al.l {2005) and IHesman et al.l (j2007f l. pro- 
duc es a two-level CO pro file that is roughly similar to 
the ILellouch et al.l (|2005T ) CO vertical profile (see Ta- 
ble 2]), although still outside of their quoted uncertain- 
ties; and a tropospheric C O abundance that is consistent 
with IHesman et al.l (]2007l ) . Using a cooler thermal pro- 



file also pushes the best-fit transition pressure to lower 
pressures (Fig. \T§§. As an additional test of the ef- 
fect of the thermal pr ofile on the CO best- fit solution, 
we used t he plot from IHesman et al.l (|2007l ) to approx- 
imate the IHesman et al.1 (|2007f ) data, and modeled the 
CO (3-2) line in a similar way to our own data. We 
note that we do not inc lude the small leakage correc- 
tion (|Hesman et al.ll2~007ft when converting from antenna 
temperature to brightness tempera ture. We find that our 
best-fit model to the approximate Hes man et al.l (|2007l ) 
data using our nominal thermal profile has a CO mole 
fraction of 1.1 ppm in the upper atmosphere; this value 
increases to 1.7 ppm of CO in the upper atmosphere if 
we perform the fit using a cooler th ermal profile, which i s 
closer to the best-fit value found bv lHesman et al.l (2007). 
This behavior is similar to what we see when fitting our 
own CO (2-1) and (1-0) line data. 

While we find better agreement with the CO p rofiles 
of ILellouch et al.l (|2005D and IHesman et al.l (|2007l) when 
we use a cooler thermal profile, our results do favor the 
warmer nominal thermal profile. As Figs. l2"0"landl2"T1 illus- 
trate, mod el spectra prod u ced u sing the ILellouch et all 
( 2005) and IHesman et"aT1 1)20071 ) CO and thermal pro- 
files, but with our radiative transfer code, are less con- 
sistent with our data than our nominal fit; this is par- 
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Fig. 20. — Comparison of the CO (2-1) line data (this work) 
with a selection of best-fit CO profiles from previous authors 
(see Table [7]). In each case, we use the reported CO and tem- 
perature profiles from the indicated study. Using our model 
code, we fit the 1- and 3- mm gain uncertainties to mini- 
mize x 2 and plot the scaled model. For reference, the best-fit 
two-level solution from this work is plotted as well (black). 



ticularly true in the line wings, 0.5-2 GHz from line cen- 
ter. The Hesman et al.1 (|2007l ) solution also appears to 
be too high in the CO (2-1) emission peak. In addi- 
tion to their use of a cooler thermal profile, the remain- 
ing cause of discrepa ncy between our result and that of 
ILellouch et all ( 2005) may be the more limited frequency 
coverage of their dataset: if we consider only our data 
within the frequency range covered by the ILellouch et al.l 
(2005) dataset, we find a best-fit vertical CO profile with 
1.00 ± 0.08 ppm in the stratosphere (above 3 mbar) and 
0.62 ± 0.05 ppm in the troposphere. We considered two 
additional factors that may contribute to the difference 
betwe en our solution a nd that found bv | Hesman et all 
(|2007f k the fact that iHesman et"aT1 (|2007f ) observe a 
higher frequency line, which probes slightly higher al- 
titudes in the atmosphere (Fig. [T]); and the uncertainty 
in the JCMT beam effi ciency, which is used to convert 
the lHesman et all ([2007) data from antenna temperature 
to brightness t empe rature. By fitting the approximate 
IHesman et al.l (|2007t ) data, we find that of the two, the 
uncertainty in the JCMT beam efficiency has a stronger 
effect on the best-fit CO profile solution. The authors al- 
low for a 5% uncertainty in the beam efficiency. With this 
constraint, our best fit to their data has a mole fraction 



Fig. 21.— Same as Fig. [20l except for the CO (1-0) line 
data. 



of 0.5 ppm in the lower stratosphere and troposphere. If 
we allow for a beam efficiency unc ertainty of ±10%, our 
best-fit models to the approximate IHesman et all (2007) 
data have mole fractions of only 0.0-0.1 ppm of CO in 
the lower atmosphere. We therefore conclude that an 
underestimate of the uncertainties in the beam efficiency 
could potentially account for the remaining difference be- 
tween our best-fit s olution and the best-fit solution of 
IHesman et ail (|2007l) . 

Two other recent papers look at the CO abundance on 
Neptune using (far)infrared data. iFletcher etldl (|2010f ) 
observed fluorescent lines with AKARI, and fit two pro- 
files to the CO (2-1) fluorescent line: in their Profile 1, 
CO is limited to altitudes a b ove 1 mbar. Profile 2 is 
similar to the He sman et al.l ([20071 1 best-fit profile: 2.1 
ppm of CO is present in the stratosphere, decreasing by 
a factor of four at altitudes below 10 mbar. They find 
that they require some CO below 10 mba r (their Profile 
2) to reproduce their data. As with the IHesman et ahl 
(120071) result, we find that the IFletcher et al l (|2010l) CO 
profile produces models that are a poor fit in the line 
core and far wings. We note that for a physical CO dis- 
tribution, the CO abundance will not be constant with 
altitude, and the derived value for the stratospheric CO 
abundance will be dependent on the pressure levels one 
is sensitive to (see Fig. [§]). Usi ng Herschel m e asure - 
ments of CO lines at 153-187 /xm, ILellouch et all pOlOh 
fin d a stratospher i c CO abundance that is similar to 
the ILellouch et all (|2005l) number; roughly 1 ppm. This 
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TABLE 7 

Goodness-of-fit x 2 values obtained using the best-fit CO profiles from previous authors with our radiative transfer code and 

our data. See Section [7] 



CO profile reference 


TP profile 


X 2 (l mm) 


X 2 (3 mm) 


X 2 all 


this work^ 


nominal 


15.3 


14.6 


15.3 


Marten e~al. (2005) b 


Marten et al. (2005) 


70.2 


71.5 


70.3 


this work^ 


nominal 


10.7 


9.3 


11.1 


Lellouch et al. (2005) d 


Bezard et al. (1991) e 


15.8 


19.9 


16.3 


Hesman et al. 2007J. 


Hesman et al. (2007) 


18.0 


21.0 


18.4 


Fletcher et al. (201uV 


nominal _ 


25.3 


16.8 


24.0 


Lellouch et al. (2010V 


Lellouctfet al. (2010) 


14.2 


16.6 


14.5 



a best-fit one-level model, from Table [3] 

1.0 ppm CO everywhere 
c best-fit two-level model, no H2S, from Table [4] 

d 1.0 ppm CO at altitudes above 20 mb ar; 0.5 ppm CO deeper than 20 mbar 
e same as used in Lellouch et al. (2005) 

*2.2 ppm CO at altitudes above 6 mbar; 0.6 ppm deeper than 6 mbar 
g 2.1 ppm CO at altitudes above 10 mbar; 0.5 ppm deeper than 10 mbar 

h from lFletcher et al.l |20TCD 

'we assume the same CO profile as in Lellouch et al. (2005); only the stratospheric CO was constrained from Lellouch et al. (2010) 



value is consistent with our re sults (though inconsistent 
with the lFletcher et al.l (|2010f ) solution). 

7.2. Implications: internal CO 

The CO abundance that originates from Neptune's 
deep atmosphere acts as a probe of Neptune's global 
oxygen abundance: according to the net thermochemi- 
cal reaction (Eq. 1) the equilibrium CO mole fraction 
is directly proportional to the equilibrium abundance of 
H2O, and under the conditions of Neptune's deep atmo- 
sphere, nearly all the gas p hase oxygen is conta i ned i n 
water. As first described by iPrinn and Barshavl (|1977t ). 
the observed tropospheric CO mole fraction (0.0-0.3 ppm 
from our analysis) represents the equilibrium abundance 
at the CO 'quench level', which is defined as the depth 
at which 

T~chem — ^mix (7) 

where r c ^ em is the timescale for chemical conversion 
of CO into CH4 and T m i x is the atmospheric mixing 
timescale. Above the quench level, vertical mixing trans- 
ports CO to higher altitudes before the constituents have 
a chance to equilibrate {T c h em > i~ m i X ); and the CO mole 
fraction remains constant at a level which can be much 
higher than the equilibrium value. The H2O mole frac- 
tion can then be determined from the CO mole fraction 
via the equilibrium relation 



g CO (gR 2 ) 3 (P T ) 2 
qCR 4 qR 2 



AjG u (CO)-A / G u (CH 4 )-AjG u (H 2 0) 



18) 



where R is the gas constant, T is the temperature and 
Pr is the total pressure in bars at the quench level, qX is 
the mole fraction of species X, and A/G°(X) is the Gibbs 
free energy of formation of species X at the temperature 
of the quench level. 

Determining the CO quench level requires knowledge of 
the limiting reaction rate for converting CO to CH4; this 
is set by the slowest reaction step in the fastest chem- 
ical pathway for CO — > CH4 conversion. The origina l 
chemical scheme analyzed by IPrinn and Barshavl (jl977l ) 



for Jupiter had as the rate-limiting step 



H 2 CO + H 2 ^CTPj + OH 



0) 



This scheme was adopted by lLodders and Fegleyl 
(1994) to estimate the O/H ratio implied by the early 
detec ti ons of ~1 ppm of CO o n Neptune dMarten et al.1 
199H iRosenqvist et al.l Il992t iGuilloteau et al.l 119931: 
Navlor et al.l 11994ft . Thev found that a 440 times so- 
lar oxygen abundance is required to produce a 1 ppm 
CO abundance; this is roughly a factor of 10 higher 
than the C/H enrichment suggested by CH4 measure- 
ments. Using up dated laboratory measure ments for the 
rat e coefficients. iGrifnth and Yelld (|1999t ) showed that 
the|P rinn and Barshavl (|197l " reaction is actually about 
4 orders of magnitude slower that originally calculated, 
and is theref ore too slow to play a ro le in CO quench- 
ing kinetics. IGriffith and Y cllc (199 §) adopted instead 
a chemical scheme first advocated bv lYung" et~aT1 dHH) 
for their analysis of Gliese 229B, which has as its rate- 
limiting step 



H + H 2 CO + M -» CH3O + M 



(10) 



where M repr esents any third body (atom or mo l ecule) . 
iBezard et all (|2002l ) also use the lYung et al.l (|1988h 
scheme for th eir analysis of CO chemistry on Jupiter - 
Mor e recently. iVisscher et al.l (120101 ) , iMoses et all (|2011f ) 
and IVisscher and Moses! (|2011l ) incorporate further up- 
dates to reaction rate coefficients, and compare the 
rates of all relevant reactions to determine the dom- 
inant kinetic mech a nism for CO — > CH4 conversion. 
IVisscher and Mosesl (|2011f ) find that the dominant path- 
way for conditions in Jupiter, cool brown dwarfs and hot 
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Jupiters is: 



H + CO + M-s-HCO + M (11a) 

H 2 +HCO -> H 2 CO + H (lib) 

H + H 2 CO + M -4 CH 2 OH + M (11c) 

H 2 + CH 2 OH^CH 3 OH + H (lid) 

CH 3 OH + M CH 3 + OH + M (lie) 

H 2 + CH 3 ^CH 4 + H (llf) 

H + OH + M->H 2 + M (llg) 



in which Eq. (|lle[) is th e rate-limiting step. (Note: 
iVisscher and Mosesl ()2011[ ) indicate a second reaction 
scheme which may be important under some conditions. 
We include this second reaction in our models as well, 
but find the rate is always more than two orders of mag- 
nitude slower than the above scheme. For simplicity, we 

omit this second reaction pathway here.) 

A s econd source of error in the lLodders and Feglevl 
( 1994) deter mination of N eptune's CO quench level was 
described by [Smith (1993), who showed that the typical 
estimate of the mixing time 



L 2 /K 



(12) 



is not correct when the pressure scale height H is used 
for the characteris tic mix i ng len gth scale L. According to 
the calculations of lSmithj (Q998), effective mixing lengths 
are typically of order 0.1 0.2H which means that as- 
suming L — H will lead to an overestimate of T m i X by 
up to two o rders of magnitude. 

Using the IVisscher and Mosesl (|2011l ) rate-limiting step 
and the ISmithl (I 1998) recipe for estimating the effective 
mixing length, we evaluate the O/H enrichment implied 
by our observed deep CO mole fractions. 

Pursuant to Eq. [FJ the equilibrium CO abundance 
depends on the H 2 , CH 4 and H 2 mole fractions, tem- 
perature, and pressure. Therefore, in order to determine 
the CO quench level implied by our data, we extend our 
model of the thermal profile and composition to the deep 
atmosphere (P of order 10 5 bar). As described in Sec- 
tion 13.11 the atmospheric C, N and S abundances are 
constrained by previous studies. We assume that in the 
deep atmosphere the C/H and S/H enrichments are ~50 
times the protosolar value, and that the CH4 and H2S 
abundances represent the total elemental abundances of 
C and S, respectively. For nitrogen, we assume that the 
NH 3 /H ratio is equal to the protosolar N/H value, as 
constrained by cm wavelength observations, noting that 
the total atmospheric N/H enrichment can be higher if 
much of the nitrogen is in N2. The O/H enrichment is the 
quantity we wish to determine from these calculations. 

For large enrichments of H2O, the approximation 
O g as/H gH 2 0/ (2 • gH 2 ), where qX is defined as the 
mole fraction of species X, does not hold, since a signifi- 
cant fraction of the hydrogen is contained in species other 
than H2 . Therefore to relate the solar oxygen enrichment 
to the various abundances we approximate 



<?H 2 



2 • qU 2 + 2 • qR 2 + 3 • (?NH 3 + 4 • <jCH 4 + 2 • qU 2 S 



(13) 



Similar equations relate the C/H, S/H and NH 3 /H en- 
richments to the mole fractions of CH4, H 2 S and NH 3 . 
Table [8] presents the deep atmospheric mole fractions for 
relevant species, for several different values of the oxygen 
enrichment. In relating the global composition of Nep- 
tune to protosolar abundances, we have assumed that 
the atmospheric abundances of C, N, O and S repre- 
sent Neptune's global C, N, O and S enrichments. In 
particular, we ignore the removal of oxygen from the at- 
mosphere by the formation of rock: for protosolar com- 
positi on gas, of orde r 20% of the oxygen is tra pped in 
rock (jLoddersI [200l IVisscher and Mosesl |20T1 . This 
means that Eq. [T2]is an underestimate of Neptune's total 
glob al O/H ratio. The pro tosolar abundances are taken 
from I Asplund et "all (|2009f ) . We note that the values for 
the protosolar abundances have been revised significantly 
over the past two decades; for example, the protosolar 
O/H ratio ha s been adjusted downward s by more than 
10% from the[G revesse and Noelsl (|1993l ) values used by 
lLodders and Feglevl (| 19941) . 

The model used for calculating the thermal profile 
down t o high pressures is described in Ide Pater et al.l 
(|1991a| ) and references within; starting with the deep 
composition (Table [8j at an initial temperature and pres- 
sure, the model calculates the temperature at successive 
steps upwards in altitude by assuming the temperature 
follows either (1) a dry adiabat (regardless of the rela- 
tive humidity of the atmospheric constituents); or (2) the 
appropriate wet adiabat (Fig. |22|) . We expect the wet 
adiabat case to be more appropriate. Clouds condense 
when the partial pressure of a given species exceeds its 
saturation vapor pressure, and the atmospheric compo- 
sition changes accordingly. No clouds form at temper- 
atures above the critical temperature of water (647 K). 
To calculate the adiabatic TP profile, we use an ideal 
gas equation of state and "intermediate" hydrogen, i.e. 
the specific heat is near that of normal hydrogen, and 
the ortho to p ara ratio is close to the equilibrium value 
(|Wallacd [r980). As in Section 13731 we adjust the start- 
ing temperature at the deep est laye r so th at the thermal 
profile roughly matches the iLind al (1992|) value of 71.5 
K at 1 bar. 

Utilizing the IVisscher and Moses (2011) rate-limiting 
reaction, the chemical lifetime of CO is given by 



7~chem 



[CO] 



[CO] 



d[CO]/dt frna[CH 3 OH][M] 



(14) 



where [X] is the number density of species X and fyiq e | 
is the reaction rate of step (|lle[) . Given the rate of the 
reverse reaction £[XTe]R we can write T c hem as 
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TABLE 8 

Mole fractions in the deep atmosphere (below condensation levels) for different O enrichments. The enrichment factor for C is 

48 x solar, S is 50 x solar, N is lx solar, and the O enrichment varies. 



Molecule 


mole fraction 
O/H=50x solar 


= 100 x solar 


=400 x solar 


=600 x 


solar 


=700 x solar 


H 2 


0.785 


0.744 


0.483 


0.292 




0.191 


He 


0.141 


0.134 


0.0869 


0.0526 




0.0345 


CH 4 


0.0250 


0.0252 


0.0266 


0.0276 




0.0281 


NH 3 


1.28 x 10" 4 


1.29 x 10" 4 


1.36 x 10" 4 


1.41 x 


10~ 4 


1.44 x 10 -4 


H 2 


0.0473 


0.0954 


0.402 


0.626 




0.744 


H 2 S 


1.28 x 10 -3 


1.29 x 10" 3 


1.36 x 10 -3 


1.41 x 


10 -3 


1.44 x 10 -3 
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Fig. 22. — Thermal profiles for the deep atmosphere, assum- 
ing an adiabatic extrapolation. The value at 1 bar is matched 
to the Voyager 2 measurement. The red profile assumes a dry 
adiabatic lapse rate throughout the atmosphere; the blue pro- 
file assumes the appropriate wet lapse rate at expected cloud 
condensation levels. The range in the profiles indicates vari- 
ation in composition, from 50 to 600 times the solar O/H 
ratio. The range of CO quench temperatures (T qU ench) for 
K = 10 7 - 10 9 cm 2 s 1 is indicated by the vertical grey bar. 



of temp erature are taken from the NIST-JANAF tables 
(jChasel[T998h . The re verse reaction rate i P h is cal- 
culated as described in lVisscher and Moses (|2Llilf l using 
the expression 



where 



0=1 + 



l + (fco[M]/fcoc) 

log 10 (fc [M]/fcoo 

0.75- 1.271og 1Q . 



(19) 



(20) 



and the para meters for calculating ^ fTTefc are given by 
I Jasper et all (|2007| ): 



k = 1.932 x i 3 T - 9 - 88 e- 7544 / T 

+ 5.109 x io-"T- 6 - 25 e- 1433 / T cm 6 s" 1 

fcoo = 1.031 X l -10 r -0 018 e 16.74/T cm 3 g -l 

F c = 0.1855e- T / 155 ' 8 + 0.8145 e - T / 1675 + e " 4531 / T 



(21) 
(22) 
(23) 

(24) 

We calculate the atmospheri c mixing time scale T m i X 
according to Eq. (fl"2f using the iSmithl (|1998l) recipe for 
determining the effective mixing length L. In general, 
we find th a t L w 0.12 — 0.18H, which is consistent with 
the ISmithl (|1998T ) findings. The eddy mixing coefficient 
in the deep atmosphere is constrain ed by the observed 
heat flux <fi = 433 ± 46 erg cm~ 2 s _1 (Pearl a nd Conrathl 
119911). Using the scaling relationship for free dry convec- 
tion |Stoni[l97|) 



7~che 



[CO] 



[CHgHOHHM]^ 
nK eq 

^x[H 2 ]2[M]to 

K eq 



(15) 
(16) 
(17) 



where Pt is the total pressure in bars, n is the total num- 
ber density in cm~ 3 , and K eq is the equilibrium constant 
of OH + CH 3 o- CO + 2H 2 : 



K eq = e 



4, G°(CO)-A , C°(CH,)-A, G°(OH) 

— 1 (18) 



The Gibbs free energies of formation as a function 



K . 



(25) 



we estimate that K is of order 2 x 10 8 cm 2 s — 1 . In Eq. 
(|25|) . p is the gas density, R is the gas constant and 
Cp is the specific heat at constant pressure. Following 
lLodders and Feglevl (|1994l) we assume a plausible range 
for K of 10 7 - 10 9 cm 2 s" 1 . 

Using Eqs. ([7|) and ((HJ, we now solve for the CO quench 
level and the predicted CO mole fraction, as a function 
of the eddy diffusion coefficient and the assumed value 
of the O/H ratio. We find that, for K = 10 7 - 10 9 cm 2 
s _1 , Eq. {7J becomes true at a temperature T quenc h = 
850 — 1100 K; coincidenta lly this is quite similar to the 
lLodders and Feglevl (|1994f) value of T quench = 998 K. The 
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quench temperature, which is indicated in Fig. I22[ is in- 
sensitive to the thermal profile (dry or wet lapse rate) 
used. Since the quench level is below the condensation 
level of the deepest cloud, the abundances of H 2 , CH4, 
and H2O are the same for the dry and wet adiabat cases; 
however the pressure at which the quench temperature 
is reached changes (Fig. |2"2"|) . which affects the CO mole 
fraction at the quench level (Eq. [8J. The quench level 
(and therefore, the predicted observable) CO mole frac- 
tions are plotted in Fig. [23] for the dry and wet adiabat 
cases; for reference, we have indicated CO abundances 
of 0.03, 0.1 and 0.3 ppm with horizontal lines. We find 
that, for the dry adiabat case, 0.1 ppm of upwelled CO 
and K < 10 9 cm 2 s _1 implies a global oxygen enrich- 
ment of at least 400 times solar; this is roughly 8 times 
the C/H enrichment implied by Neptune's observed CH4 
abundance. For a wet adiabat, 0.1 ppm of upwelled CO 
and K < 10 9 cm 2 s^ 1 implies a global oxygen enrichment 
of at least 650 times solar. 

We might expect the planets to be uniformly enriched 
in heavy elements (jOwen "eFaLl[T999h : if this is so, then 
for fast mixing (K~ 10 9 cm 2 s _1 ) we would expect an 
upwelled CO abundance of 1.5 x 10 -9 or 8.4 x 10~ n , 
for the dry and wet adiabat cases, respectively. We note 
that such low internal CO mole fractions are consistent 
with our data. 

7.3. Implications: external CO 

In addition to the abundance of CO that is vertically 
transported from Neptune's deep atmosphere, a signif- 
icant external source (4>co = 0.5 — 20 x 10 s molecules 
cm~ 2 s _1 ) is implied by our analysis (see Table[5]). These 
inferred CO production rates are in gen eral agreement 
with t he 10 s molecules cm~ 2 s _1 found bv lLellouch et all 
(2005) and lower than the 10 10 molecules cm -2 s _1 es- 
timated by iHesman et al.l (|2007) . The precise value of 
4>co depends strongly on the eddy diffusion and thermal 
profiles used in the model. 

The high rates of 0co inferred by the stratospheric 
CO mole ratio, and the high CO/H2O ratio on Neptune 
have led to the suggestion that comet impacts are the 
primary supp ly mechanism of CO to Neptune's uppe r 
stratosphere (jLellouch et al.l 120051: IHesman et al.l 120071) . 
We explore effectiveness of producing these abundances 
of CO through cometary impacts by calculating the ap- 
proximate CO production from esti mates of impac t rates , 
similar to the analysis performed by B ezard et al.1 ()2002[ ) 
for Jupiter. We use for the impact ra t es of comets at 
Jupiter the estimate from iZahnle et all pOOl : 

Nj(D > 1.5 km) = 0.005lg;oo3 y r_1 ( 26 ) 

as well as their approximation that the rate for Neptune 
is roughly 



N N (D > 1.5 km) » 0.25iV, 7 (D > 1.5 km) 



(27) 



where D is the comet diameter. We test two possible 
comet size distributions: the first is the Shoemaker and 
Wolfe (1982) distribution: 



A>jv(> D) = N N (D > 1.5 km) 



D(km) 
1.5 



Secondly, we test the IZahnle et al.l d2003) "Case B" im- 
pact distribution, which was determined from the size 
distribution of craters on Triton: 



iVjv(> D) = 2.62 N N {D > 1.5km) ( — -j— ^ ) yr" 1 



1.5 

(D < 1.5 km) (29a) 



N N {> D) = 0.129 N N {D > 1.5km) [ - ■ ) 



-2.5 



1.5km) ( K ' 

(D > 1.5 km) (29b) 



Using these comet distributions, we calculate the mass 
rate of CO due to comet impacts as 

M = /o/co / 4nR 2 pN(> R) dR g yr" 1 (30) 

JRmin 

where R is the comet radius, /o is the fraction of 
the comet mass in the form of oxygen and fco is the 
fraction of the oxy gen that ends up as CO. Following 
IBezard et al.l (|2002f ) we adopt values of /o = 0.5 and 
/co — 0.9, and a typical comet nucleus density p = 0.55 
g cm~ 3 . 

T o calculate th e size range of comets to use, we fol- 
low IBezard et al.1 ((2002) in choosing a minimum radius 
of H m in = 0.15 km, and calculate the maximum radius 
Rmax as the maximum size for a comet for which the 
impact timescale is equal to the diffusion timescale: 



N{> D max ) 



TK 



The diffusion timescale can be approximated as 



tk 



2H$ 
K 



(31) 



(32) 



yr- 1 (28) 



where Ho and Ko are the scale height and diffusion rate 
at the tropopause, where mixing is slowest. We consider 
three values of K : 200, 800 and 2000 cm 2 s" 1 . The 
slowest mixing rate of 200 cm 2 s _1 corresponds to the 
minimum eddy diffusion rate in the Romani 'A' profile 
and our 'test' eddy profile; the Romani 'B' and Moses 
'2' profiles have Ko values near 2000 cm 2 s _1 , and the 
Moses '1' profile has an intermediate Ko ~ 800 cm 2 s _1 . 
The CO production rate 4>co is determined from the CO 
mass influx rate M for the range of Nn(D > 1.5 km) and 
for the two different comet size distributions described 
above. These calculated CO production rates are shown 
in Table M 

Using a CO flux to the upper atmosphere of 0.5 — 
20 x 10 s cm~ 2 s -1 , as inferred by our physical fits, 
we calculate the diameter fli of a single large comet 
impact required to reproduce the observed CO abun- 
dance, given tk and for each value of Ko. We also 
report the timescales for impacts of size D\ - which 
can be compared to the t imescale estimates re ported by 
iLellouch et ail (j2005ft and lHesman et al.l (|2007l ). We note 
that, without detailed knowledge of the vertical eddy 
diffusion profil e , the two-l evel CO profil e s re ported by 
ILellouch et~"aLl (|2005l ) and IHesman et"ail (|2007l ) do not 
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Fig. 23. — Predicted CO mole fraction in the visible atmosphere of Neptune due to upwelling from the deep atmosphere as a 
function of the eddy mixing coefficient K and the O/H enrichment over solar, using a dry adiabatic extrapolation (left) and a 
wet adiabatic extrapolation (right) for the thermal profile. Horizontal lines indicating 0.3, 0.1 and 0.03 ppm of CO are shown 
(red); the shaded region indicates the plausible range of K values. O/H enrichments of 50, 100, 400 and 600 times solar are 
shown. 



preclude a constant CO influx by (sub)kilometer-sized 
comets (rather than the single event discussed by these 
authors.) Eddy diffusion profiles like the Romani 'B' and 
Moses '2' profiles will produce CO profiles that, when 
mimicked by a two-level profile, will have transition pres- 
sures in the 10 mbar range (Section 16.41) . Since very 
large impacts are rare, a constant infall rate of smaller 
comets may be a preferable explanati on. We find that, o f 
our two comet size distributions, the lZahnle et al.l (|2003f ) 
Case B distribution gives more optimistic predictions for 
CO production: the rates of CO injection found for the 
small and intermediate Ko values are ^3 x 10 s 

cuT 2 s^ 1 and O^Oj^;^ x 10 8 cm~ 2 s _1 , respectively, 
suggesting that cometary impacts are in fact a feasible 
mechanism for supplying the observed abu ndance of CO 
to Nep tune's upper atmosphere. However, IZahnle et al.l 
(2003) suggest that their impact rate for 1.5 km comets 
at Neptune is likely an overestimate, meaning the true 
CO production rate due to comets could be smaller. 

8. SUMMARY/CONCLUSIONS 

We have observed Neptune in the CO (2-1) and (1-0) 
rotational lines with the Combined Array for Research 
in Millimeter-wave Astronomy. Our r adiative trans- 
fer an alysis indicates a preference for the iFletcher et al.l 
(2010) thermal profile in Neptune's upper atmosphere 
over warmer and cooler profiles, and we find that the 
best-fit solution for the CO vertical profile is strongly 
dependent on t he atm ospheric TP profile. Adopting the 
IFletcher et al.l (|2010D thermal profile, we find that good 
fits to the data (both to the individual lines and to the 
combined dataset) are characterized by CO mole frac- 



tions of O.IIq' 2 parts per million (ppm) in the tropo- 
sphere, and l.llg'g ppm in the stratosphere. Higher CO 
mole fractions, particularly in the stratosphere, are fa- 
vored when cooler thermal profiles are used. If the CO 
mole fraction resulting from vertical mixing is greater 
than 0.1 ppm, our calculations imply an O/H abundance 
that is at least 400 times the protosolar value. How- 
ever, since we do not rule out a 0.0 ppm tropospheric CO 
mole fraction, we cannot place a lower limit on Neptune's 
global O/H ratio, and our data do not require that Nep- 
tune's oxygen enrichment exceeds its carbon enrichment. 
We find that the stratospheric deposition rate of CO is 
0.5 — 20 x 10 s CO molecules cm~ 2 s _1 ; such a high abun- 
dance could be supplied by impacts from (sub)kilometer- 
sized comets, as long as the eddy diffusion rate near the 
tropopause is small. 

This work will be followed up with spatially-resolved 
mapping of Neptune in the CO (2-1) line with CARMA, 
to look for latitudinal variations in the CO abundance. 
ALMA will play a key role in the three dimensional map- 
ping of trace species on Neptune, as well as on the other 
Solar System giant planets. Such studies will provide 
further insight into the composition, atmospheric circu- 
lation, and environments of the giant planets. 

The data presented in this work were obtained with 
CARMA. Support for CARMA construction was derived 
from the states of California, Illinois, and Maryland, 
the James S. McDonnell Foundation, the Gordon and 
Betty Moore Foundation, the Kenneth T. and Eileen L. 
Norris Foundation, the University of Chicago, the As- 
sociates of the California Institute of Technology, and 
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TABLE 9 

Production rate of CO from comets. 











Shoemaker 


and Wolfe 




Zahnle et al. 


Case B 




Kq (cm 2 s 1 ) 


tk (yr)t 


Di 


(km)^ 


Dmax (km) 


log0_^( cm -2 s^ 1 ) 




D max (km) 


log c (cm~ 2 s -1 ) 


i£>, (yr) 


200 
800 
2000 


930 
230 
93 


3.4 
2.2 
1.6 


- 12 

- 7.4 

- 5.4 


0-81S1 
0.5±g;l 


7 i +o"5 

' • i -o.7 

6.3±°:2 


1900 - 120000 
750 - 49000 
410 - 26000 


r, , + 0.9 

z '°-0.7 
i q+0.6 

1>d -o.5 

7+ ' 4 
u -'-0.3 


7 7+U.t) 
'•'—0.6 
7 ,+0.6 
'-°-0.8 


1100 - 130000 
350 - 41000 
160 - 19000 



a diffusion time scale 

b diamctcr of single comet impact required to produce <j>CO = 0-5 — 20 X 10 s 
c total CO flux expected based on impact rates 
d l/N N (> D 1 ) 



the National Science Foundation. Ongoing CARMA 
development and operations are supported by the Na- 
tional Science Foundation under a cooperative agree- 
ment, and by the CARMA partner universities. This 
work was supported by NASA Headquarters under the 
NASA Earth and Space Science Fellowship program - 



Grant NNX10AT17H; and by NSF Grant AST-0908575. 
The authors would like to thank L. Fletcher for providing 
his temperature and CH4 profiles, and G. Orton for pro- 
viding his H2 CIA absorption coefficients. The authors 
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